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Abstract. In equilibrium systems, the conservation of the number of particles (or 
, mass) leads to the equalization of the chemical potential throughout the system. Using 

a non-equilibrium generalization of the notion of chemical potential, we investigate the 
influence of disorder and of the balance of mass fluxes on the generalized chemical 
potential in the framework of stochastic mass transport models. We focus specifically 
on the issue of local mesurements of the chemical potential. We find that while local 
dynamical disorder does not affect the measurement process, the presence of large-scale 
, geometrical heterogeneities (branching geometry) leads to unequal local measurement 

results in different points of the system. We interpret these results in terms of mass 
flux balance, and argue that the conditions for the global definition of the chemical 
k*" | potential still hold, but that local measurements fail to capture the global theoretical 

^| ■ value. 

in 

O ! 1. Introduction 

The question of defining relevant macroscopic control parameters in non-equilibrium 
systems still remains largely open [HE]- At equilibrium, thermodynamics provides us 
with intensive parameters like temperature, pressure and chemical potential, that are 
uniform throughout the system even in the presence of heterogeneities, and that can 
in most cases be easily measured. In non-equilibrium steady states, it is natural to 
try to find similar types of parameters, and different generalizations of equilibrium 
concepts have been proposed either through statistical approaches often related to 
entropy notions [3j HJ El [7], or using generalized fluctuation-dissipation relations in 
theoretical PEl[T0l[ni[T2l[l3l[Tl[T5]and experimental [HJ [13 [HJ H9] contexts - A few 
studies [HI EDI EU E2] have more explicitly considered the question of the equalization of 
such parameters throughout inhomogeneous systems, but this issue has not been settled 
yet. Such a question has also been addressed in a previous work [23j EH; where the 
equilibrium definitions of intensive thermodynamic parameters conjugated to conserved 
quantities have been extended to some classes of non-equilibrium steady-state systems. 
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Provided that a condition called "asymptotic factorization property" holds [23 [23], this 
approach yields for the out-of-equilibrium chemical potential A the definition^ 




where Z(M) is a generalized partition function. Roughly speaking, the asymptotic 
factorization condition is satisfied when only short range correlations are present in the 
system. It holds in particular in the absence of correlations, when the joint probability 
distribution is factorized (apart from the global conservation law). 

The main properties of this generalized notion of chemical potential have been 
studied in [23], where it was shown to exhibit some interesting equalization properties 
in some classes of inhomogeneous non-equilibrium systems. Possible difficulties arising 
from the non-equilibrium nature of the systems considered and especially from the 
dynamics at the contact between different systems have also been outlined. 

In the present work, we study whether the globally- defined chemical potential given 
in Eq. ([I]) can be evaluated through local measurements, and whether measurements 
performed on the system at different locations yield the same value, which is a non-trivial 
issue when the system is inhomogeneous. We also investigate how the flux balance -a 
strong constraint specific to non-equilibrium situations- may influence the results of 
local measurements of the chemical potential. As simple examples of inhomogeneous 
non-equilibrium systems, we consider different models belonging to the class of mass 
transport models introduced in Refs. [26j EZ]> i n which a globally conserved mass is 
transferred between neighbouring sites. We specifically consider both the case of local 
disorder in which the dynamics is locally heterogeneous (Section [2]), and the case of a 
large-scale geometrical heterogeneity related to the presence of several branches in the 
system (Section [3]). 

2. Mass transport model on a ring with local disorder 

2.1. Definition of the model 

We focus here on a one-dimensional mass transport model with periodic boundary 
conditions, that is on a ring geometry. On each site % = 1, . . . , N resides a positive 
mass ?Ti j. The continuous time stochastic dynamics, which preserves the total mass 
M = Y^=i m i-i is defined as follows (see Fig. []]). An amount of mass p, is transferred 
from site i, containing the mass m;, to site % + 1 with a probability per unit time 
p<Pi(n\mi), and to site i — 1 with a probability per unit time q<pi(p,\mi), where q — 1 — p 
(note that N = and N + 1 = 1 due to periodic boundary conditions). The rate 
(Pi(fi\mi) is defined as 



| Note that this generalized definition of the chemical potential differs by a factor — 1/T (where T is 
the temperature) from the conventional equilibrium definition [25] . 
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where v(/j,) and fi(m) are arbitrary positive functions, with fi(m) possibly site- 
dependent. Thus transport is biased (except if p — q), which generates a flux of mass 
along the ring. With the above rate, the steady-state distribution takes the form [26] 

1 n / n \ 



mi 



M 



(3) 



with S(x) the Dirac delta function, and where the partition function Z(M) is defined as 
Z{M) = [ f[[d mi /(mi)] 5 (j2m - m) . (4) 



" i=l \i=l / 

Note that the function v(/S) does not influence the steady-state distribution, but only 
the dynamics. When performing numerical simulations, we use throughout the paper 
v(fi) = 1. 
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Figure 1. Schematic drawing of the dynamics of the mass transport model. A fraction 
fi of the mass m, situated on site i is transferred according to the local conditional 
rates tpi(fj,\mi), either to site i + 1 with probability p, or to site i — 1 with probability 
q = 1 - p. 



2.2. Globally- defined chemical potential. 

Let us here consider for /i(mj) the simple form fi(rrii) = m* _1 with rji > for all i. To 
calculate the generalized chemical potential corresponding to the conserved mass in the 
system we need to find the dependence of the partition function Z on M: 



r N / N \ 

Z(M)= Yl [dnii m?" 1 ] 6 ^ m, - M 
•* i=i \i=i j 



(5) 



where the integrals are over the positive real axis. A simple rescaling = x^M reveals 
the searched dependence: 

jv r N ( N \ 

z{m) = me^i*- 1 / ni^i^ -1 ] m - 1 

J i=l \i=l / 

= D N M N *>- 1 (6) 

with rj = A r_1 Y^iLiVii an d where is a constant independent of M. The generalized 
chemical potential is obtained from the derivative of In Z 
d\nZ Nrj-1 

X ~ ~dM ~ ~W~ (7) 
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leading in the thermodynamic limit N — > oo to 

A=5, (8) 
where p = M/N denotes the average density. 



2. 3. Local measurement of the chemical potential. 

Once the chemical potential has been theoretically defined, an important issue is to 
know whether it can be measured. One possible way to perform a measurement is to 
connect to the considered system a probe system, with a much smaller size in order 
not to significantly perturb the main system. The probe system is assumed to have a 
known equation of state, so that its chemical potential can be deduced from its average 
mass. The connection between the two systems is in general local -just like when a 
thermometer is put into contact with a system to measure its temperature. As a result, 
it is important to verify that the measurement result is the same wherever the connection 
is made. This issue becomes non-trivial when the system is inhomogeneous. 




Figure 2. Sketch of the contact of a small probe system (right) with a large system 
(left). Mass is transferred in the direction of the arrows with the indicated transport 
rates (see text). 



To check this feature, we consider the above inhomogeneous mass transport model 
on a ring with transport rates <p;(/z|mj) defined according to Eq. fl2]) with fi(m) = m Vi ~ l . 
We choose a sinusoidal space dependence for iji, of the form rji = 2 + sin(27n/iV), 
% = 1,...,N. We also set p = 1 in the numerical simulations. As a probe system, 
we use a homogeneous mass transport model, with site-independent rates wM 777 '*) 
defined with / pr (m) = 1, that we successively attach to sites i = N/4 and i = 3iV/4 of 
the inhomogeneous system (see Fig. [2]). 

Fig. E] shows that the chemical potential X(t) of the probe converges to the chemical 
potential of the main system (computed from its equation of state). The results are seen 
to be independent of the location where the probe is attached to the inhomogeneous 
system. We further observe that even in a completely disordered system, in which 
the rjiS are independent and identically distributed quenched random variables, the 
measured value of the chemical potential is still the same as long as the mean value 
rj — l/NJ2iVi remains the same as in the deterministic case. 
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Figure 3. (Color online) Chemical potentials A(t) of the probe system (upper curves) 
and X(t) of the main system (lower curve) during the measurement process, in three 
different situations: (i) Inhomogeneous system with rji = 2 + sin(27ri/iV), and probe 
connected to i — N/A (red curve); (ii) Same system, with the probe connected to 
i = 3-/V/4 (green curve); (iii) Inhomogeneous system with random and uncorrelated 
values of r\i on each site, drawn from a uniform distribution on the interval 1 < rji < 3 
(blue curve). In all cases, the mean value rj = 2 is the same. The contact is switched 
on at time t = (vertical dotted line). Initial densities arc p = 2 for the main system 
and p pr = 0.2 for the probe. All measurements converge to the theoretically expected 
value (dashed line). System size TV = 15870, probe size N pr = 512. 




2-4- Relation with the balance of mass flux. 

We have seen that, in this simple model with a ring geometry, the locally-measured 
chemical potential coincides with the global one, in spite of the heterogeneities which 
makes the local density site-dependent. This result can be interpreted from the one-site 
probability distribution p^rrii). Considering the rest of the system as a reservoir of 
mass, one finds for the probability distribution on site i 

Mm l ) = ±-f l (m)e- Xm % (9) 

where Qi is a normalization constant. Hence the local chemical potential is indeed the 
same everywhere in the system. 

Interestingly, this result can also be given an alternative interpretation in terms of 
flux balance. Starting from Eq. (Q, one can relate the flux and the chemical potential 
A. The total mass flux $j crossing site i can be expressed as 

poo rm 

®i = (p-q) dmpi(m) dnfitpi(ti\m). (10) 
Jo Jo 

which can be rewritten, using Eqs. (EJ) and (j2J), as 

$i = (p-q) dnnvMe-* (11) 
Jo 

One thus obtains that $j = $ is independent of the site considered, as expected from 
the steady-state flux balance in this linear geometry. Alternatively, one could interpret 
Eq. ( 1TT]) by saying that the uniformity of the local chemical potential results from the 
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flux balance. Hence one can guess that flux balance plays an important role in the 
determination of the local chemical potential. This role will appear even more clearly 
in the example considered in the next section, where the geometry is no longer purely 
linear. 

3. Model with three branches 

3.1. Definition of the model 

We now investigate the issue of the local measurement of the chemical potential in a 
model where the local dynamics is essentially homogeneous, but where the heterogeneity 
results from a branching geometry. We consider a mass transport model with three 
branches, corresponding to the geometry displayed in Fig. HI The three branches are 
assumed to be oriented. The transfer rate from site % to the neighbouring site is p ip(fj,\rrii) 
along the positive direction (according to the orientation of each branch) and qip(fj,\mi) 
along the negative direction. The rate y?(/x|m) is defined according to Eq. (j2j), with 
a site-independent weight function f(m). At the branching points (A — > B or C), 
probability rates for the transfer to branches B and C are reweighted by factors jb and 
7c, as shown on Fig.Hl For instance, at the lower branching point on Fig. HI the transfer 
from branch A occurs with rate 7spy 9 (/ i l m i) to branch B, and with rate 7cP I m i) t° 
branch C. Due to this specific geometry, the probability distribution does not necessarily 
factorize, even with the choice of transport rates ip(fi\m) given in Eq. (j2J). 
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Figure 4. Sketch of the model with three branches. Left panel: Mass is transported 
along each oriented branch according to the rate p (p(n, m) in the positive direction, and 
to the rate q <p(fi, m) in the negative direction. Right panel: Zoom on the branching 
points, where specific rules are taken into account, some of the transfer rates being 
reweighted by factors jb and 7c- 



For mass transport models defined on an arbitrary graph, a sufficient condition for 
the factorization of the probability distribution has been given in Ref . [28] . In the most 
general case, such models are defined by transport rates tpij(p,\mi) from site i to site j 
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of the form 



him) 



(12) 



where the function Vij(fi) is identically zero if there is no directed link from % to j. If 
the condition 



holds for every site z, the probability distribution factorizes, with the local probability 
weight given by /i(mj) [28], as in Eq. (J5J). 

In the case of the model with three branches, the branching sites violate the 
sufficient condition for factorization given Eq. ( fl3l) . which is a strong indication that 
strict factorization does not hold. Hence the exact solution of this model is not known, 
and it is not clear a priori whether the present model satisfies or not the "asymptotic 
factorization condition" , which is a key criterion for the existence of a globally defined 
chemical potential. Indeed, although the strict factorization property is likely to be 
violated, one might ask whether this violation is 'localized' around the branching points 
and if the local chemical potential in the bulk of the branches remains uniform. It is thus 
interesting in this situation to perform numerical measurements of the local chemical 
potential. 

3.2. Measure of the chemical potential with a probe system 

To test this issue, we measure the chemical potential with a probe system attached 
to a bulk site of the branch considered, in analogy to the procedure explained in 
Sect. 12.31 The transport rates are the same in all branches as well as in the probe 
system, tp(^\m) = 1. All three branches have the same number of sites JV&, so that the 
system size is iV = 3A^ (to be specific, the two branching points are included in branch 
A). The branching points contribute to the number of sites in branch A. Simulations 
are done with p — 1 and Jb — Jc = h- The results of the numerical implementation 
of the measurement are shown in Fig. [5J The locally measured chemical potential is 
denoted as \ u in branch v = A, B, C . It turns out that we obtain equal values for the 
locally measured chemical potential in branches B and C, but a very different value for 
branch A. 

3. 3. Measure of the chemical potential through local fluctuations 

To check whether this result obtained by probing the system depends on the 
measurement method, we apply in the following an alternative measurement technique 
that does not require any external device. The idea, introduced in an earlier work |24j . 
is based on a direct measurement of the fluctuations in small subsystems of different 
sizes, within each branch. We briefly sketch the procedure in the following. Let us first 
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Figure 5. Chemical potentials A„, v = A, B, or C measured with a probe system 
attached to branches A, B and C, plotted as a function of time t after the contact is 
switched on at t = 0. The transport rates are the same in all branches and in the probe 
system, ip(fx\rrii) = 1. The dynamics at the contact is defined as shown in Fig. [2] The 
initial value of the global density is p = 0.5 in the main system and p pr = 0.2 in the 
probe system. All three branches have the same size Nf, = 43348, and the probe size 
is iVp r = 1024. Measurements on each branch have been performed during separate 
runs. 



define the quantity g v (v = A, B, or C) as the variance of the total mass M v in branch 
v divided by the number of sites N u in this branch: 



We have checked numerically that g v does not depend on the size of the subsystem 
chosen within a given branch. Consequently g v is an intensive quantity depending only 
on the local density, g v = g v {p v ). For a detailed description of this procedure see [21]. 

We choose the same parameters for the dynamics on the three branches as in 
Sec. 13.21 Each branch is of the same size Nf, and the transport rates are given by 
(p(fj,\m) = 1 everywhere. The result of a numerical implementation of the measurement 
of the function g v (p v ) in the three-branch model is shown in the left panel of Fig. [6j 
The functional behaviour of g v {pv) is seen to be the same for the three branches, and we 
denote this function simply as g{p u )- We further observe on Fig. [6] that the numerically 
measured g(p u ) is very close to the theoretical value g(p u ) = p 2 corresponding to a 
homogeneous system with the same local dynamics ip(p\m) = 1, so that we shall use 
this theoretical expression in the following. 

Using the grandcanonical ensemble [21], the chemical potential can be related to 
the function g(p u ) according to 




(15) 



from which the relation 




(16) 
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follows. This means that the functional behaviour of the local chemical potential A with 
density p„, that is the local equation of state, is the same in the three branches, namely 
A(p^) = p" 1 for the specific dynamics chosen here. This result is not surprising since the 
local dynamics is the same everywhere and although the probability distribution does 
not factorize, only weak correlations are expected in the bulk of each branch. Taking 
into account the fact that the local densities p„ in the three branches are not equal, it 
follows that the values X u = A(p„) of the local chemical potential of the three branches 
differ as well: A a ^ \b = Ac- 





^~'o6o6666556556 



0.2 



O.ft 



O.N 



Figure 6. Left panel: measure of the function g v {p v ) in each branch A, B and C. 
The same functional form, very close to the theoretical form g(p) = p 2 corresponding 
to a homogeneous system, is obtained in each case. Transport rates: <p(p\m) = 1. 
All branches have the same size Nf, = 341. Right panel: locally- measured chemical 
potential A, y = A(p v ) in each branch v — A, B, and C, plotted as a function of the 
total density in the system p = M/N (from simulations). The expected value of A„ for 
ip(p,\m) = 1 is X v = p^ 1 (see text). 



Though this discrepancy points to a weakness of the concept of non-equilibrium 
chemical potential, it can however be understood by taking into account the balance of 
fluxes. Starting from the expression (II ip of the flux we first note that if the function 
v(jj) in the transport rates is identically 1, v(fi) = 1, then the flux is given by $ = 1/A 2 
(note that the expression of the flux is independent of the form of f(m)). We use here the 
local chemical potential A, as the flux is governed by the local dynamics -see Eq. ( TTTj) . 
We shall come back to this point in Sec. 13.41 In the bulk of the branches, where the 
correlations are expected to be small, the different fluxes are given by $^ = 1/A 2 A and 
&c — &b — l/A^ (fluxes in branches B and C are equal since 7b = 7c)- Moreover 
we know that the flux in branch A is twice the flux in branch B, which leads to the 
following relation between the chemical potentials of the different branches: 

Aa = ^ = ^- (17) 

This relation has been verified numerically for tp(/j,\m) = 1, where we expect analytically 
\ v = p^ 1 , v — A, B,C - see the right panel of Fig. |6] and left panel of Fig. [7] for the 
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numerical results. We can generalize this flux balance argument to situations where 
branches B and C do not have the same flux, in which case the relation reads: 

x * = Xb = J¥?« ■ (18) 

Note that at equilibrium, when the fluxes vanish, this correction is not present and all 
chemical potentials equalize. But as soon as a little bias is introduced in the dynamics, 
the value of the chemical potential in branch A differs from that of the two other 
branches. Let us emphasize that this difference is not perturbative with respect to the 
bias: — A# does not go to zero when p — q — > 0, but rather remains constant as long 
as p > q. 




U U,t U,U U,0 1 ~ 5 ' 1 * " 

P "P 

Figure 7. Left panel: rescaled values A^ r ^, denned as = A^ and X^ c = \b,c/V%, 
plotted as a function of p, showing a collapse of the data (same data as in the right 
panel of Fig. |6]) . Right panel: differences in the locally- measured chemical potential 
of system B and A (x), and system C and A (+) as a function of p compared to the 
analytical prediction a(~p) = 1/(3/5 (black line) as discussed in Sect. 13.51 

From a more theoretical perspective, the discrepancy between A^, A# and \ c 
questions the validity of the asymptotic factorization condition [231 121] required to 
define a global chemical potential from Eq. ([Q). It would thus be interesting to verify 
explicitly whether this asymptotic factorization property holds or not. This is a difficult 
task because (to our knowlegde) the exact joint probability distribution is not known for 
generic rates in the present geometry, as it does not fulfill Eq. ( 1T3|) . Yet, the specific case 
of the Zero Range Process (ZRP), where masses are discrete, turns out to be solvable 
and thus deserves to be investigated in more details. 

3.4. A solvable case with discrete masses 

The ZRP case, which has been intensively studied in the literature [29], corresponds to 
the choice v (fi) = 5(fi — 1), so that masses irij take integer values, denoted as Hi in the 
following. We study here the ZRP in the three branch geometry illustrated on Fig. HJ 
The dynamics is defined as in Sec. 13. 1[ except that can only take the value /1 = 1. To 
simplify the calculations, we set 7b = 7c = 1 and choose p > q. 
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The relation fllip between flux and local chemical potential reads in this case 
$„ = e~ Xv . Given that fluxes are different in the three branches, it is clear that 
the chemical potentials are not equal. More quantitatively, the balance of fluxes 
$j4 = + $<7 implies Xa = — ln($^/$ B ) + A#. From this expression we see that 
there is now a shift, given by the logarithm of the ratio of the fluxes, in the value of the 
chemical potential. 

To better understand the origin of this shift, we compute the steady-state 
distribution of the ZRP with three branches. Interestingly, it has been shown that 
the steady-state distribution of the ZRP on an arbitrary graph remains factorized [28J . 
Slightly rephrasing the results of Ref. [28], one finds that the distribution P({ni}) is 
given by 

where 5 is the Kronecker delta symbol, and where the local "fugacities" Zi satisfy for all 
i the equation 

H V ij Z i=Yl V Ji Z 3 ( 20 ) 

with = %(1) defined in Eq. f|T2|) . In the present three-branch model, is equal 
to p if there is a positively oriented link from % to j, and to q if the link is negatively 
oriented; = in the absence of link. Note that the fugacities Zj are defined only up 
to an overall arbitrary factor. In the present model, condition ( l20l) reads, for all site i 
different from a branching point: 

pZi-i -Zi + qz i+1 = 0. (21) 

The general solution of this linear equation (valid separately on each branch) is a linear 
combination of solutions of the form Zj = r J , with a parameter r obeying the relation 

qr 2 -r + p = 0. (22) 

Eq. (1221) has two solutions, T\ — 1 and r 2 = p/q > 1. As a result, Zj can be expressed 
on each branch as a linear combination of the two independent solutions rj and r J 2 , 

N v -f 

(23) 



Z 3 




where ~z v and K v are constants [y = A, B, or C) and N v denotes the number of sites in 
branch v. The total number of sites in the system is given by N = Na + Nb + Nc + 2 
(the last term accounts for the two nodes). Taking into account Eq. ( 1201) . formulated for 
the two branching points, one can match the expressions (123]) corresponding to different 
branches and determine the constants K v as well as the ratios za/~zb and za/zc- In 
the limit where the sizes iV^, Nb and Nc of the three branches go to infinity, one 
finds za/z b = za/zq = 2. As Zj is defined up to an overall prefactor, one can choose 
for instance 1b = ~zc — 1 an d Za = 2. Given that the terms proportional to K v are 
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exponentially decaying corrections, it turns out that Zj is almost constant, and equal to 
z v , on each branch v. 

Hence the definition of the chemical potential proposed in [23j El] can indeed be 
applied. Splitting the system into two parts, branch A on one side, and branches B and 
C on the other side, one obtains 

dlnZ A d\nZ BC 
OM A oM BC 
with Ma = J2ieA n i an d M B c = J2ieB,c n i t ne respective masses of the two subsystems. 
Since Za — 2, we get, neglecting the exponential corrections appearing in Eq. (123]) 

Z A = 2 M -Z A °\ Z BC = zf c , (25) 

where Z A ^ and Zg C are the standard partition function of these two subsystems, that 
would be obtained by taking Zi = 1 for all i. This results, in the large size limit, in 

A A = ln2 + A A , \ BC = \ BC . (26) 

We thus recover in this way the result directly obtained from the balance of fluxes, 
showing that this result is actually consistent with the asymptotic factorization condition 
on which the definition of the chemical potential is built (let us recall that the probability 
distribution of the present ZRP model is factorized). Note that the shift In 2 is 
independent of the precise value of the bias p — q, as long as this bias is non-zero. If 
p = q, equilibrium is recovered, and the bias vanishes. The presence of this bias is thus 
a genuine non-equilibrium effect, which appears non-perturbatively, in a discontinuous 
way. 

It is also interesting to note that the shift In 2 cancels out from the one-site mass 
distribution. Indeed, one has for a site i in the bulk of branch A, using Eq. (1260 



Piin) = ^- fin) 2 n e~ XAn = ±- f(n) e~^ n (27) 

and for a site in branch B or C 

Pi(n) = -^—f(n)e-~ XBcn . (28) 

Hence both distributions have exactly the same form as a function of the local chemical 
potential, making the shift undetectable from a local measurement. 

3. 5. Consistency between flux balance and equality of the chemical potentials 

We have seen in the ZRP case that the discrepancy between the local chemical potentials 
X u can be explained by the presence of exponential factors that differ from one subsystem 
to the other. It is tempting to try to generalize this scenario to the case of the mass 
transport model with continuous masses, although no analytical solution is available in 
this case. At a heuristic level, a more general form of Eq. fl25|) can be proposed, namely 

Z A = erW M *Z$\ Z BC = Z { B d l (29) 
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since the global average density p is fixed, and can thus a priori enter as a parameter in 
the exponential factor. Eq. (126]) is then replaced by 

Aa = «(p) + X A , X B c = X B c- (30) 

Assuming that the asymptotic factorization condition holds, one has A^ = Xbc, so that 

a(p) + X A = X B c- (31) 

The constant a(~p) is determined as follows. Given a value of p, one looks for the densities 
Pa(p) and psc(p) satisfying the contraints 

$( P a) = 2$(pbc), (32) 
p = Kp A + (1 - k)pbc (33) 

with k = Na/{Na + Nbc)-, an d where $(p) is the value of the local flux for a local 
density p. The parameter a(~p) is then obtained, consistently with Eq. ( I3TT) . as 

a(p)=A(p B cr(p))-A(pA(p)). (34) 

In the ZRP case, one recovers a = In 2, independently of the density p, while with 
continuous masses, assuming v(p) = 1 and f(m) = m r, ~ l , one finds 

3hi — l\ . , 

(35) 

This result is in agreement with the numerical findings of Sect. 13.31 where rj = 1 and 
k — |, yielding a = l/(3p) (see right panel of FigJTJ). Hence, despite the absence of 
an exact solution of the mass transport model (beyond the specific ZRP case), this 
tentative scenario provides a consistent explanation of the numerical results obtained 
for the local chemical potential, showing that the asymptotic factorization condition 
may still hold while local measurements fail to find an equilibrated chemical potential 
throughout the system. 




4. Discussion and conclusion 



In this paper, we have attempted to clarify the influence of the flux balance condition on 
the generalized chemical potential in non-equilibrium mass transport models. We have 
seen in particular that even when the asymptotic factorization condition holds, locally 
measured values of the chemical potential may not equalize, while the theoretically- 
defined ones remain equal. This surprising property has been traced back to the presence 
of exponential factors appearing in the probability weights of different subsystems. 
These exponential factors cancel out at the level of the local statistics (for instance 
the single-site distribution) and are thus locally undetectable. But on the other hand 
these factors play an essential role in the global statistics to ensure the balance of flux. 

Quite surpringly, such factors can exhibit a discontinuity between equilibrium and 
weakly non-equilibrium situations, as exemplified by the exact solution of the ZRP: at 
equilibrium (zero flux), no exponential factor is present (and local chemical potentials 
equalize), while the tiniest flux generates a factor 2 M , leading to a shift In 2 between 
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the local chemical potentials. We have also shown that a tentative generalization of 
this scenario, assuming a density- dependent shift a(p), is consistent with the numerical 
results. 

An important consequence of these results is that in most cases it is not possible 
to define along the present lines a locally measurable chemical potential that would 
equalize throughout the system. This comes from the fact that the shift a(p) depends 
on the global density of the system, which cannot be measured locally. Hence, although 
no long-range correlations are present, a global information is required to relate what 
happens at two different locations in the system. 
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